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ABSTRACT 

We consider the effect of radiation pressure from ionizing photons on black hole (BH) mass estimates based 
on the application of the virial theorem to broad emission lines in AGN spectra. BH masses based only on the 
virial product AV 2 R and neglecting the effect of radiation pressure can be severely underestimated especially 
in objects close to the Eddington limit. We provide an empirical calibration of the correction for radiation 
pressure and we show that it is consistent with a simple physical model in which BLR clouds are optically 
thick to ionizing radiation and have average column densities of Nh ~ 10 23 cm" 2 . This value is remarkably 
similar to what is required in standard BLR photoionization models to explain observed spectra. With the 
inclusion of radiation pressure the discrepancy between virial BH masses based on single epoch spectra and on 
reverberation mapping data drops from 0.4 to 0.2 dex rms. The use of single epoch observations as surrogates of 
reverberation mapping campaigns can thus provide more accurate BH masses than previously thought. Finally, 
we show that Narrow Line Seyfert 1 (NLS 1) galaxies have apparently low BH masses because they are radiating 
close to their Eddington limit. After the radiation pressure correction, NLS 1 galaxies have BH masses similar 
to other broad line AGNs and follow the same MBH-Ce/Lsph relations as other active and normal galaxies. 
Radiation forces arising from ionizing photon momentum deposition constitute an important physical effect 
which must be taken into account when computing virial BH masses. 

Subject headings: radiation mechanisms: general — galaxies: active — galaxies: fundamental parameters — 
galaxies: nuclei — quasars: emission lines — galaxies: Seyfert 



1. INTRODUCTION 

In the last few years, it has become increasingly clear that 
supermassive black holes (BH) are an essential element in 
the evolution of galaxies. The key observational evidence 
of a link between a BH and its host galaxy is provided 
by the tight correlations between BH mass and luminosity, 
mass, velocity dispersion and surface brightness profile of the 
host s p heroids (pCormendv & Richstondll995llGebhardt et al.1 
2000; iFerrarese & Merritj l2000t iMarconi & Hunt) 120031 
Graham & Driverl 120071) . The link between BH and host 



galaxy is probably established by the feedback of the ac- 
creting BH, i.e. t he active galactic n ucleus, on the host 
galaxy itself (e . g. ISilk & Reesl 11998b iGranato et all l2004t 
IDi Matteo etafl 120051; ICroton et all 120061 and references 
therein). 

In order to fully understand the implications of BH growth 
on the evolution of the host galaxies it is fundamental to 
measure BH masses in large samples of galaxies from zero 
to high redshifts. Direct BH mass estimates based on stel- 
lar and gas kinematics are possible only in the local uni- 
verse and their comp lexity does not allow their application 
to lar ge samples (e.g. IFerrarese & Fordll2005l; iMarconi et all 
120061) . The limit of the local universe can be overcome 
with the reverberation mapping (RM) technique (see, e.g., 
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iPeterson & Bent z 2006 for a recent review) which provides 
an estimate of the Broad Line Region (BLR) average dis- 
tance from the BH C^blr)- The BH mass can thus be de- 
rived using the virial theorem, Mbh = /AV 2 Rblr/G, where 
AV is the width of the broad emission line and / is a scal- 
ing factor which depends on the p hysical properties of the 
BLR (e.g. lPeterson & Wandell2000l) . Although this technique 
is potentially plagued by man y unknown systematic errors 
dKrohldl200lUCollin et all 120061) . BH masses from reverber- 
ation mappping are in a greement with the MsH -Ce relation 
of normal galaxies (e.g. lMcLure & Dunlopll2002i) . However, 
this technique is very demanding in terms of telescope time 
and it can be applied only to a few objects espe cially at high 
redshifts (Peterson et al. 2004; Kaspi et aT1l2007l). The radius- 
luminosity relation discovered by Kas pi et all ([2000) shows 
that continuum luminos i ty can be used as a proxy for /?blr 
(iKaspi et al.ll2000t 12005b iBentz et al.ll2006bl) . From the spec- 
trum of a broad line AGN it is therefore possible to obtain a 
single epoch (SE) BH mass estimate. 

One of the most important sources of uncertainty i n viria l 
Mbh estimates is the scaling factor /. lOnken et al.1 (120041) 
have provided an estima te of / assum i ng th at the AGN 
in the RM database of IPeterson et~all (120041) follow the 
M B H-o" e relation of no rmal galaxies dTremaine et al.1 120021: 
IFerrarese & Fordll2005l) . The factor / by lOnken et all (120041) 
is only app licable to estimates of the virial product based 
on RM (see Peterson et al. 20 04] for more details). Building 
on the results by Onken et al. (2004), Vestergaard & Peterson 
(2006) have calibrated scaling relations for SE virial Mbh es- 
timates which combine the width of broad H/3 with the lumi- 
nosities of XL\ at 5100A. 

Overall, SE virial estimates are commonly used to estimate 
BH masses in large sample of galaxies from zero to high 
redshifts (e.g. IWillott et al.1 120031 iMcLure & Dunlopl 120041; 
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Vestergaard 200H lJiang et~aT1l2007l) and are deemed accurate 
only from a statistical point of view on large samples of ob- 
jects since a sin gle measurement can be wron g even by a fac- 
tor of ~ 10 (e.g. lVestergaard & Petersonll2006l) . 

There are three important considerations which are sug- 
gested by the results presented in the above papers. First, SE 
virial BH masses of a few objects (e.g. high z, high L quasars 
or Narrow Line Seyfert 1 galaxies) imply they radiate near or 
above the Eddington limit. The virial theorem is based on the 
assumption that the system is gravitationally bound and this 
might be violated in super-Eddington sources where the out- 
ward force due to radiation pressure overcomes gravitational 
attraction. Second, even when L < L^i, one should take into 
account that the radiation force partially compensates gravi- 
tational attraction. In the standard accretion disk model, the 
source of ionizing photons can be c onsidered point-like at the 
distance of the BLR (see however Col lin & Hurell200ll for a 
different point of view) and the radiation force scales as r~ 2 
mirroring the radial dependence of the BH gravitational at- 
traction. Thus BLR clouds are effectively being pulled by 
a smaller effective BH mass and all present virial mass es- 
timates for objects close to their Eddington limit, where ra- 
diation pressure is not considered, might be underestimated. 
Finally, the Eddington limit is computed assuming that the ra- 
diation pressure is due only to Thomson scattering of photons 
by free electrons. As supported by reverberation mapping, 
by the radius -luminosity relation an d other observational ev- 
idences (e.g. iBlandford et al.| [r990). BLR clouds are almost 
certainly photoionized. Thus BLR clouds are subject to ra- 
diation forces arising from the deposition of momentum by 
ionizing photons which can substantially exceed that due to 
scattering. 

The importance of radiation pressure due to ionizing pho- 
tons and its possible effects on virial BH masses has already 
been mentioned in a few papers (e.g. lMathewil99l iGaskell 
1996) but seem to have not been considered in detail subse- 
quently. This effect might be particularly important in Nar- 
row Line Seyfert 1 galaxies which are believed to accrete 
close to their Eddington limit. Indeed, they are character- 
ized by small BH masses com pared to other A GNs and to the 
MBH-£sph/o" e relations (e.g. iMathur et al.ll2001l) . It has also 
been noted that the distance of NLS1 galaxies from the Mbh- 
L S piJcr e relations is larger for objects with larger Eddington 
ratios dGrupe & Mafhuri 12004) suggesting that smaller BHs 
are growing faster. Alternatively, this might be an indication 
that that virial BH mass are underestimated in the high L/Lndd 
regime. 

In this paper we investigate the effect of radiation pressure 
on virial BH mass estimates. In §|2]we present a simple phys- 
ical model for the radiation pressure effect on virial BH mass 
estimates. In §[3] we calibrate the effect of radiat ion pressure 
on vir ial B H masses adapting the p r ocedu res of lOnken et al.l 
d2004l) and lVestergaard & Peterson! (120061) . In §g]we apply 
our corrected virial BH mass estimates to Narrow Line Seyfert 
1 galaxies and show that these galaxies are indeed consistent 
with the MBH-Ce/^sph relations, showing that BHs are not ab- 
normally small. Finally, we discuss our results and draw our 
conclusions in §|5] 

2. THE EFFECT OF RADIATION PRESSURE ON VIRIAL BLACK 
HOLE MASS ESTIMATES: A SIMPLE PHYSICAL APPROACH 

We will explore the effect of radiation pressure on BLR 
clouds using a simplified model which assumes that (i) each 
cloud is optically thick to ionizing photons but optically thin 
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FIG. 1 . — The ratio between virial BH masses taking into account radiation 
pressure (Mbh) and those based only on the virial product (Mbh.o) is plotted 
against the classical Eddington ratio based on Mbh.o ■ Mbh/^bh.o is thus the 
correction factor which should be applied to BH mass estimates based only 
on the virial product. A?h is, on average, the total column density of each BLR 
cloud along the direction to the ionizing source. 



to scattering processes, (ii) the Thomson cross-section is rep- 
resentative of all scattering processes involving free or bound 
electrons and (iii) both recombination and scattered photons 
are 'isotropically' re-emitted. These assumptions are valid 
if Uc/a B (H) < Nh < l/er T where Nn is the total cloud col- 
umn density along the direction to the ionizing source, U is 
the ionization parameter, chb(H) is the 'case-B' recombination 
coefficient for hydrogen and or is the Thomson cross-section. 
For typical condi tions in the BLR (T e ~ 2 x 10 4 K, U ~ 0.01, 
e.g. lNetzerll2006l) 1.2 x 10 21 cm" 2 <JV H < 1.5 x 10 24 cm" 2 . 

The total force acting on a cloud in the outward radial di- 
rection and due to radiation pressure is 

r+oo r 

F = l d^(l-e-)AA (1) 

where L v is the luminosity of the AGN continuum emis- 
sion, r is the cloud distance from the ionizing source, t v is 
the optical depth of absorption/scattering processes and AA 
is the cloud surface exposed to the AGN radiation. Scattering 
is important only for non-ionizing photons therefore, under 
the above assumptions, it is possible to write: 

~ Lion A , L '-'ion , T * . 

F=— r AA+ a T N H AA (2) 

where the two terms are the radiation forces due to absorp- 
tion of ionizing photons and Thomson scattering, respectively, 
Lion is the total luminosi ty of the AGN ionizing con tinuum, 
hv > 13.6eV (see, e.g.. lPetersonlll997l:lKroliklll999l) andjV H 
is, on average, the total column density of each BLR cloud 
along the direction to the ionizing source. The contribution 
to the radiation force from the absorption of line photons is 
negligible for the opticall y thick clouds cons idered here (see, 
e.g., the seminal paper bv lCastor et alj|1975l) . 
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Taking into account the total radiation force acting on each 
cloud and assuming that the BLR is a bound system, it is pos- 
sible to derive a modified version of the classical virial the- 
orem which takes into account radiation as well as gravita- 
tional forces. Approximating the cloud mass as ~ m p A^AA 
the modified expression for the virial BH mass Mbh is: 



BH : 



V z r L 
f-pr + J 



\-a+- 



M, 



(3) 



where / is a geometrical factor which takes into account the 
geometry of the BLR, LEdd.o is the classical Eddington lumi- 
nosity for a solar mass object and a = L; on /L. This expres- 
sion has a physical meaning as long as the system is bound, 
i.e. as long as the radiation force on BLR clouds is smaller 
than gravity. In formulae 

^Edd.O 



L< 



[l-a + a/(a T N H )] 



(4) 



where LEdd.o is the classical Eddington luminosity. Neglect- 
ing momentum injection by ionizing photons (a = 0) we re- 
cover the classical relation L < Lndd.o- Using Mbh from 
Eq.[3]to compute LEdd.o, it should be noticed that for L — > oo, 
L/LEdd,o — ■> 1/ [l-a + a /(ctjNh)] and L/L£dd,o will always be 
less than or equal to 1 . This is a consequence of the assump- 
tion of gravitationally bound BLR which allowed us to write 
Eq. [3] Therefore, it is not possible to establish whether a sys- 
tem is above Eddington by using virial BH mass estimates, 
since they are themselves based on the assumption of a sub- 
Eddington system. 

In order to quantify the effect of the radiation force correc- 
tion we write Eq.[3]as 



M BH = M BH ,o 



1 + - 



-Edd.O 



l-a + - 



(5) 



where Mbh.o is the standard virial BH mass computed with- 
out taking into account radiation pressure. In Fig.Q]we show 
the behavior of Mbh/Mbh.o as a function of L/LEdd,o and for 
different values of A^h- The a = Lj 0n /L bolometric correction 
has been computed following Marconi et al. ( 2004) and is on 
average a ~ 0.6 in the 10 10 — 10 I2 L Q luminosity range. For 
A^H=10 23 cm~ 2 andL/LEdd,o > 0.1, Mbh/Mbh.o varies between 
2 and 10. This can be much larger for smaller column densi- 
ties of BLR clouds but values at low A^h should be taken with 
caution since the adopted formula is valid only if the cloud 
is optically thick to ionizing photons, i.e. A^h > Uc/as(H) ~ 
1.2 x 10 2I (£^/0.01)cm~ 2 . The correction factor remains small 
(< 2) only for column densities A^h > 10 24 cm" 2 . Clearly the 
correcting factor critically depends on the A^h value which sets 
the total cloud mass and thus the relative importance of grav- 
itational attraction with respect to radiation pressure. Over- 
all, this figure suggest that neglecting the effect of radiation 
pressure might result in Mbh values which are underestimated 
even by a factor ~ 10. 

Virial estimates of BH masses are based on the assump- 
tion that the BLR is gravitationally bound to the BH and that 
outflowing motions are negligible. In recent years, building 
upon observational evidence for outflows in the BLR, alter- 
native models have been proposed i n which part of the BLR 
is in the form of a disk wind (e.g. iMurrav & Chiangl 1 1 995b 



Chiang & Murray! 11996b Elvis 2000b ICollin & Hurd l200lb 
Proga et ai1l2000b iProgd 12007b lEverettl 120051 and references 
therein). This possibility has gener ated a debate about the re- 
liability of virial BH masses (e.g. iPeterson & Wandelll2000b 



i Krolikl I200H lOnken & Petersonl 12002b ICollin et al.1 12006b 
Vestereaard & Peterso nl l2006l) which is beyond the scope of 
this paper. Nevertheless, virial BH mass estimators are widely 
used, and in order to investigate the effect of radiation pres- 
sure on such estimates, we must necessarily start from the 
same set of assumptions for our simple model. 

3. THE EFFECT OF RADIATION PRESSURE ON VIRIAL BLACK 
HOLE MASS ESTIMATES: AN OBSERVATIONAL APPROACH 

The simple physical approach presented in the previous sec- 
tions suggests that virial BH mass estimates can be written as 
a function of observed quantities as 



; VlR ^ ( L 5100 

G \ 10 44 ergs 1 



M, ; 



(6) 



where L5100 represents \L\ at 5100A. After Eq. [3] g corre- 
sponds to 



g = 6.0x 10 6 



9.0 



\-a + - 



(JtNh 



(7) 



where b = L/L5100 is the bolo metric correction at 5100A. Fol- 
lowing iMarccHfietal] (120041) . the L/L5100 bolometric correc- 
tion is on average b ~ 9.0 in the 10 10 - 10 12 Lq luminosity 
range. / and g are free unknown parameters which depend on 
the physical and geometrical properties of the BLR. In partic- 
ular the g factor critically depends on the assumed A^h value 
which determines the cloud mass and thus sets the relative 
importance of gravity and radiation pressure. 

A correction for the radiation force which is proportional 
to L is more general than the simple physical model pre- 
sented in the previous section therefore, in order to avoid 
a-priori assumptions on the values of the physical param- 
eters characterizing BLR clouds, we can determine / and 
g following a procedur e simil ar to lOnken et alj (120041) and 
IVestergaard & Petersonl (120061) . Thus our model will only 
provide a simple physical interpretation of the empirical g val- 
ues. 

3.1. Black hole masses from reverberation mapping data 

Onk en et all (|2004) considere d the AGNs f r om th e re- 
verberation mapping database by IPeterson et al.l d2004l) with 
measured stellar velocity dispersion. They used the time lag 
of the broad lines for R and the velocity dispersion of the 
r.m.s. spectra for V. They determined / by assuming that the 
AGNs in their sample follow the MBH-a" e relation for normal 
galaxies. 

We first update the RM database by IPeterson et alJ (120041) 
with the newer esti mates of BLR time lags for NGC 4151 
dBentz et al l l2006al). NGC 4593 dDennev et alj 120061) and 
NGC5548 dBentz et al.ll2007h . We exclude from the database 
PG121 1 +143 and IC4329A b ecause their time lags are not 
reliable (IPeterson et al.l 12004). When possible, we correct 
the average AGN luminosities ALa(5100 A ) for th e host 
galaxy contamination following iBentz et all d2006bl). We 
consid er the host galaxy velocity dispersions by Onk en et alJ 
(120041) and we supplem ent them with the da ta by lNelson et alJ 
d200l for Mrk 279 and lDasyra et"ai1 J2007h for PG1229+204, 
PG1426+015, PG1617+175 and PG2130+099. 

/ and g are then derived by finding the minimum of: 



X 



[(logM B H),-(logMBH)0,/] 2 

(5 logM BH ) 2 + (6 lo S M BH ) 2 0i + AS 2 



(8) 
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TABLE 1 
Calibration of RM virial masses 



TABLE 2 
Calibration of SE virial masses 



Database 



A rc; 



Onken2004t 
Onken2004 
Updated 
Updated (Faml) 
(Fam2) 



5.5 (+1.9;-1.5) 

5.2 (+1.6;-1.2) 
4.8 (+1.5;-1.3) 
3.1 (+1.3; -1.5) 

4.3 (+1.2;— 1.1) 



-10* 
-10* 
-10* 
7.6 (+0.3;- 

< 2 



■0.3) 



0.39 
0.39 
0.52 
0.50 



t AS = 0.0 as in Onken et al. 1 2004).* Fixed fit parameter. 

where (logMBH); is the log BH mass of the z'-th object which 
depends on / and g, (logM B H)o,i = a + /31og(cr e /200kms _1 ),- 
is the expected mass value from the M ^-a relation 
dTremaine et alj|2002t iFerrarese & Ford 2005]). a e is the stel- 
lar velocity dispersion of the host spheroid, (5\ogM^)i is 
the error on (logMBH); based on the errors on V 2 , R and 
^(logg e ); is the error on log(<7 e /200),;. At variance with 



Onk en et al.l ((2004), we allow for an intrinsic dispersion of 
the Mbh-c relation, AS, which we assume equal to 0.25 dex 
(e.g jTremaine et al l 12001 iMarconi et all 120041: iTundo et al l 
2007). We follow a standard x 2 minimization and we es- 
timate errors on the para meters with the bootstrap method 
dEfron & Tibshiranill 19941) with 1000 realizati ons of the par- 
ent sample. As shown b vl Onken et all ( 2004) the use of the 
IFerrarese & Ford! d2005l) or lTremaine et al.l d2002l) version of 
the MBH-c e relation provides consistent results; therefore, in 
the following we will focus only on the lTremaine et al.l (120021) 
relation, a = 8.13 ± 0.06, /3 = 4.02 ± 0.32. 

The results of the fitting procedure ar e summarized in Ta- 
ble Q] We have considered the original lOnken et al.l (120041) 
database and the updated one. Errors on fit parameters are 
determined from the percentiles of the bootstrap results at the 
68% confidence level around the median. Several consider- 
ations can be made from the results in Table Q] As a sanity 
check, we are able to reproduce the results by lOnken et al.1 
(2004) i.e. / = 5.5 ± 1.9 (first table row). The fits shown in 
the second and third row indicate that when g is fixed and 
negligible, the use of the updated database or the use of an in- 
trinsic dispersion for Mbh-c do not significantly change the / 
value. With the use of the updated database which has a larger 
number of objects, the scatter of the residuals is significantly 
increased. When g is free to vary, the bootstrap analysis shows 
that there are two distinct families of solutions: those where 
both / and g are determined and those where g is negligible 
and totally undetermined. The existence of two families of 
solutions from the bootstrap simulations is an indication that 
the dependence on luminosity can be inferred only from part 
of the sample, i.e. from the objects with the largest L/Lem ra- 
tios. In roughly 20% of the sample realizations the number 
of these objects is low, g is undetermin ed and the / values are 
consistent with the Onke n~eTai1 d2004l) determination. The in- 
clusion of the g parameter has the net effect of decreasing /, 
since the expected BH mass is fixed by the M^-a^ relation. 

Our ability to determine an accurate empirical value of g 
is limited, as were previous efforts to determine /, by the 
size, composition and accuracy of the existing reverberation 
database. In particular, it currently contains few sources with 
high Eddington ratios, which provide the tightest constraints 
on g. With this caveat in mind, however, we provide a first es- 
timate of / = 3.1 ± 1.4 and logg = 7.6 ±0.3 to compute Mbh 
from reverberation mapping data. 



Mbh from RM 




log/ logg 




A res 




M BH from SE (V H/3 , L 5l00 ) 


/ = 5.5;logg = 
/ = 5.5;logg = 
/ = 3.1;logg = 


-iot 

-10 
7.6 


6.93 (+0.12;-0.13) -10* 
6.47 (+0.17;-0.22) 7.48 (+0.16; 
6.13(+0.15;-0.30) 7.72 (+0.06; 


-0.25) 
-0.05) 


0.43 
0.34 
0.22 


M BH from SE (V H /3, L HI3 ) 


/ = 5.5;logg = 
/ = 3.1;logg = 


-10+ 
7.6 


6.69 (+0.12;-0.08) -10* 
5.95 (+0.12; -0.17) 7.82 (+0.07; 


-0.09) 


0.46 
0.27 



t Original Onken et al. 12004) database.* Fixed fit parameter. 

3.2. Black Hole masses from single epoch spectra 

IVeste rgaard & Petersonl d2006h considered the AGNs in the 
iPeterson et alj d2004l) database. They collected single epoch 
spectra for the same sources and used the FWHM of the broad 
lines as an estimate of V and the continuum or broad line lu- 
minosity to estimat e R from the radius-luminosity relation of 
iBentz et all d2006bl) . Then they determined the corresponding 
/ parameter (see Eq.|9]below) by rescaling the virial products 
from si ngle epoch spectra to the BH masses determined fol- 
lowing |Onken!etij] ( 120041) . 

We consider the database of single epoch measure- 
ments of FWHM(H0) (herea fter V H /3), L H/3 and L5100 by 
Vestergaard & Peterson (2006) and following those authors 
we write the virial BH mass from single epoch measurements 



M BH /M =/ 



~c ( V ' 



H6 



+g 



V lOOOkms" 1 
L5 100 



( L 5 



100 



V 10 44 



ergs" 



10 44 



ergs 



(9) 



where the proxy for V is now the FWHM of the H/3 line and 
the BLR radius R is given by the radius-luminosity relation 
with a slope of 0.50 ±0.06. As before, the best / and g values 
follow from \ 2 minimization as in Eq. [8] where (logM B H)o.; 
is now the BH mass from reve rberation mapping computed 
according to lOnken et alJl2004l (/ = 5.5, logg = -10.0) or to 
our new calibration (/ = 3.1, logg = 7.6). Obviously, the AE 
term has been removed. 

The fit results are shown in Table [2] where, as before, we 
provide bootstrap errors. The fit results in the first row are 
the sanity check to show that we are able to reproduce the 
results by Vestergaa rd & Peter son (2006), who find log/ = 
6.91 ±0.02 with an rms of 0.43. Our errors are larger be- 
cause of b ootstrap simulations, but they w ould be similar to 
the ones by Vestergaard & Petersonl d2006l) if we used the for- 
mal errors of the fit. In the second row we start from the 
ass umption that virial masses from RM are computed follow- 
ing lOnken et al.l d2004l) . but we allow for a free g factor. The 
SE data are clearly able to provide an estimate of the g factor 
which turns out to be remarkably similar to what was found 
for the RM data. In the third row we start from virial RM 
masses computed with the best / and g values and there are 
two surprising results: first, the g value which turns out for 
SE virial masses is logg = 7.72 ±0.05, perfectly consistent 
with that from RM virial masses, but with a much smaller un- 
certainty. Second, the dispersion of the residuals drops from 
~ 0.4 to 0.2 dex. The latter result indicates that half of the 
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FIG. 2. — Comparison between SE and corresponding RM virial masses. Left panel: without taking into account radiation pressure (same as Fig. 8, right panel, 
of Vestereaard & Peterson 2006); the dispersion of the data along the x axis around the value is 0.4 dex. Right panel: same as right panel but taking into account 
radiation pressure as described in the text (RM virial masses are also computed with the correction for radiation pressure); the dispersion is 0.2 dex. Error bars 
combine errors on RM and SE virial masses but are dominated by the former. Error bars are different in the left and right panels because of the different relative 
importance of virial products and luminosities in RM virial masses. 

The two g values determined by (1) minimising the RM 
virial mass against the Mbh - cr e relation and (2) minimising 
the SE virial mass against the 'calibrated' RM mass are both 
consistent with a value \ogg ~ 7.7. Considering Eq. [7] we 
can derive the average Nu which is needed to obtain the g 
value determined empirically. With logg = 7.7 and a = 0.6 
we can derive Nu — 1.1 x 10 23 cm" 2 . This Nu value which 
we inferred by calibrating RM and SE virial BH masses is 
remarkably similar with the indications from photoionization 
modeling studies of the BLR. Within the framework of the 
standard BLR model, photoionization calculations can ex- 
plain observed spectra only if BLR clouds are optically thick 
to ionizing radiati on and adopted Nh are usually of the order 
of 10 23 cnT 2 (e.g. Baldwin et a Hfl99l iKasp i & Netzeri [l999t 
Korista & Goad 2004 and references therein). 



scatter of SE virial BH masses around RM ones is consistent 
with a need to take into account radiation pressure. The re- 
duced scatter of the SE virial masses is also shown in Fig. [2] 
(right panel) and should be compared with the left panel in the 
same f igure and Fig. 8 (right panel) o flVestergaard & Peterson! 

dlxjg). 

IWu et all (12004 and lGreene & Hoi d2005bl) have shown that 
it is also possible to use the luminosity of the broad H/3 in- 
stead of L5100 to avoid possible contamination of the AGN 
con tinuum emission from t he ho st galaxy. Thus, follow- 
ing ^sterga^d^&Petersoo (|2006), we substitute L5100 with 
Lh/3 in Eq. [9] to obtain SE virial masses from the broad H0 
line onl y. Inverting the L5100 relation by [Greene "&Hol 
(2005b) we can write 



Mbh/M =/ 



1000 kms" 



+g 0.732 



Lh/3 



10 42 ergs" 

0.883 



0.44 



10 42 ergs" 



(10) 



The fit results are shown in Table|2] As before, we can repro- 
duce the IVestergaard & Peterson! d2006j) calibration, log/ = 
6.67 ±0.03, and the best fit which takes into account radia- 
tion pressure shows a significant drop in the dispersion of the 
residuals providing a best fit g value which is consistent with 
previous results. 

3.3. The average column density of BLR clouds 

The results in the previous sections show that it is possi- 
ble to determine / and g both for RM and SE virial masses 
although it is difficult to accurately quantify their magnitude 
with the pr esent data. The / v alue s are smaller than those 
derive d by lOnken et all (120041) and IVestergaard & Peterson! 
(2006) because the final BH masses are still calibrated with 
the Mbh-c relation but part of the final Mbh value is accounted 
for by the effect of radiation pressure. Considering the effect 
of radiation pressure can significantly improve the agreement 
of SE and RM virial masses. 



4. THE CASE OF NARROW LINE SEYFERT 1 GALAXIES 

The nature of Narrow Line Seyfert 1 galaxies and their re- 
lation to 'normal' Seyfert 1 galaxies is still debated, but it is 
more or less generally believed that they are AGNs character- 
ized by high accretion rates and small BH masses account - 
ing for their smaller line widths (e.g. iPounds et al.l 119951) . 
Many different authors have undertaken the task of measur- 
ing virial BH masses in NLS1 galaxies and found that they 
are small com pared to br oad line AGNs with similar lumi- 
nosities (e.g. Grupe 2004). The location of NLS1 on the 
MBH-c e /L sp h, however, is still hotly debated. Most authors 
suggest that NLS 1 ga laxies have small BHs compared to their 
host galaxies (e. g. Mathu r et al.ll2001l;lGrupe & Mathurll2004t 
IZhou et al1l2006HRvan et aLl l2007) while others find an over- 
all a greement with the Msn-Q'e relation of n ormal galaxies 



a gi 

(e.g. lBotteetaill200H [Komossa & Xu 2007). A picture is 
now emerging in which the BHs in NLS1 galaxies are now 
experiencing a rapid growth which will eventually lead them 
on the MBH-^sph/oe relations as other active and normal gala x- 
ies (e.g. lCoUin & Kawaguchill2004t iMathur & Gru pe 2005). 

NLS 1 galaxies are thus ideal targets to explore the effects 
of the newly calibrated expressions which take into account 
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FIG. 3. — Distributions of BH masses for narrow line Seyfert 1 galaxies (Vu/s < 2000 kms -1 , thick line) and 'no rmal' Seyfert 1 galaxies (V H fj > 2000 km s~' , 
thin line with shaded area). Left panel: Mbh obtained with the scaling relations by Vestergaard & Peterson 12006). Right panel: A/gn obtained with the scaling 
relations which take into account radiation pressure. The numbers in the top left corners of both panels denote the mean and standard deviation of the mean 
(a/VN) of narrow and broad Seyfert 1 galaxies. 



radiation pressure. In particular, using our new calibrated ex- 
pressions for virial BH masses, we will verify (i) whether BH 
masses of NLS 1 galaxies are indeed small compared to other 
AGNs with similar luminosities and (ii) whether they lie be- 
low the MBH-tr e ,i S ph relations. 

We first test whether BH masses in NLS1 galaxies are 
on average smaller than those in 'normal' Seyfert 1 galax- 
ies. We consider the co mplete, soft X-ray selected sam- 
ple by iGrupe et al.l (2004) which is composed of 110 broad 
line AGNs with measured Vnp, Lhjs and L5100 and we com- 
pute virial BH masses using Eq. [9] In Fig. [3] we plot the 
distributions of Mbh obtained with the scaling relations by 
IVestergaard & Peterson! ((2006) (left panel) and with the scal- 
ing relations which take into account radiation pressure (right 
panel). The sample has been divided in two parts, narrow 
line Seyfert 1 galaxies (V///3 < 2000kms _I , thick line) and 
'normal' Seyfert 1 galaxies (Vh0 > 2000 kms" 1 , thin line with 
shaded area). In the top left corners of both panels we report 
the mean and standard deviation of the mean {a/ wN) of nar- 
row and broad Seyfert 1 galaxies. If radiation pressure is not 
taken into account, we recover the well known result that BH 
masses are a factor ~ 10 smaller in NLS1 galaxies. However, 
this difference is greatly reduced to a factor ~ 2 when radi- 
ation pressure is taken into account. The average BH mass 
of 'normal' Seyfert 1 galaxies is unchanged as expected since 
these objects are accreting at moderately low Eddington ra- 
tios compared to NLS 1 . It is beyond the scope of this paper to 
accurately determine the average BH mass of NLS 1 with re- 
spect to Seyfert 1 galaxies, we only wish to point out that the 
effect of radiation pressure is very important and, when taken 
into account, BH masses of NLS 1 galaxies are, on average, a 
factor 5 larger. 

We now test whether NLS 1 galaxies indeed lie below the 
^BH-c e /i sp h relations. We consider only samples where a e 
or L sp h are measured directly because we want to avoid is- 



sues connected with u sing a e surrogates like the dispersion 
of the [Om] line (e.g. iGreene & Holl2"b05at iKomossa & Xul 
1 20071). We thus c onsid er the samples of NLS 1 galaxies by 
iBotte et all d2005l) and IZhou et aL (2006) wher e a e ar e di- 
rectly measured and the sample by iRvan et al.l d2007l) . the 
only one for which accurate high resolutio n J and K photom- 
etry of the host spheroid is available. From lZhou et al.l (120061) 
we take the sub-sample of 33 sources with z <0.1 for which 
either the host galaxy appears to be face on or the SDSS fiber 
aperture is dominated by galactic bulge contribution. This 
choice is motivated by the need to avoid bulge velocity dis- 
persion values which are artificially increased by rotation of 
the galact ic disks. 

For the lBotteetaild2005h and lRvan et"aT] d2007h samples, 
we compute virial BH masses using the scaling relations by 
Vestergaard & Peterson (2006) and Eq. [9] Instead for the 
Zhou e t al. (2006) sample we use Eq. [T0l i.e. we use the lu- 
minosity of the broad H/3 as a proxy for /?blr, since, due 
the latter selection criteria, XL\ might be strongly contami- 
nated by stellar light. The comparison with exp ected BH mass 
values from the M^-a^ dTremaine et al.ll2002l) and MBH-£sph 
dMarconi & Huntll2003l) are plotted i n the Fig. B in the left 
panel w e use the virial BH masses by Vestergaard & Peterson! 

(2006) while in the right panel we use our new virial mass 
estimates which take into account radiation pressure. A more 
refined statistical analysis would be complicated by the het- 
erogeneity of the data and is beyond the scope of this paper 
but it is clear that, although with a large scatter, NLS 1 with 
'old' virial BH masses are lying preferentially below the Mbh- 
<r e relation defined by normal galaxies. When radiation pres- 
sure is taken into account in virial BH mass estimates, this 
tendency disappears or is strongly reduced. It is significant 
that th e NLS1 galaxies with bulge luminosities bv lRvan et alj 

(2007) are all lying below the expected Mbh -^sph values while 
they are in good agreement with it after radiation pressure has 
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FIG. 4. — Comparison between virial BH masses (Mbh.vu) an d those expected from the Mbh-""c or MeH-^sph correlatio ns for normal galaxies (M bh.coit. 
Tremaine et al. 2002; Marconi & Hunt 2003). Left panel: virial BH masses are computed using the calibrated relations by Vestersaard & Peterson (2006). Right 
panel: virial BH masses are computed using the relations derived in this paper which take into account radiation pressure. 



been taken into account. This is confirmed by Fig. [5] where 
we plot the histogram of the distances fr om the MeH-Ce cor- 
relation for the data by lZhou et al.l (120061) . In the top left cor- 
ner we report the mean and standard deviation of the mean 
(a/y/N) of residuals from the M^w-a e correlation. If radia- 
tion pressure is not taken into account, NLS 1 galaxies lie, on 
average, a factor ~ 5 below the correlation. However, after 
taking into account radiation pressure, virial BH masses are 
dispersed around the correlation. 

The above findings do not constitute the definitive proof 
that radiation pressure provides a solution to the small BH 
mass problem in NLS1. We only show that our calibrated 
correction for radiation pressure is approximately of the right 
amount to bring NLS1 to lie on the MBH-c e ,£ S ph correlations. 

Finally, although it is not possible to establish whether a 
system is emitting above Eddington using virial BH masses 
(see § |2), the average incre ase of BH masses by 0.5-0.7 dex 
in NL S1 galaxies (from the lGrupe et al.ll2004l and lZhou etail 
2006 samples, respectively) implies a similar decrease of their 
classical L/LEdd,o ratios. 

5. SUMMARY AND CONCLUSIONS 

In this paper we have considered the effect of radiation pres- 
sure on virial BH mass estimates. 

With a simple physical model, we have provided a correc- 
tion for the effect of radiation pressure on virial products. This 
correction mainly depends on the average column density Nh 
of broad line clouds. 

We have re-calibrated virial BH masses based on re- 
verberation mapping data and si ngle epoch s pectra fol- 
lowing a procedure analogou s to lOnken et al.l (120041) and 
IVestergaard & Peterson! (120061) . With the caveat that it is dif- 
ficult to accurately quantify the importance of radiation pres- 
sure with the present data, we find consistent values for the 
radiation pressure correction which, based on the above phys- 
ical model, indicates an average A^h ~ 10 23 cm~ 2 for BLR 
clouds. This value is remarkably consistent with the BLR 
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FIG . 5 . — Histogram of the M B H, v ii/M B H.coiT ratio for the Zh ou efafl 420061) 
sample of 33 NLS1 galaxies (s ee text). Virial BH ma s ses ar e computed us- 
ing the calibrated relations by Vestergaard & Peterson 12006) (thin line with 
shaded area) and with the relations derived in this paper which take into ac- 
count radiation pressure. 

cloud column density required in photoionization models to 
explain the observed spectra. 

When taking into account radiation pressure, the average 
rms scatter of the ratio between single epoch and reverber- 
ation mapping virial BH masses drops from 0.4 to 0.2 dex. 
The use of single epoch observations as surrogates of expen- 
sive reverberation mapping campaigns can thus provide more 
accurate virial BH masses than previously thought. 

We have considered our newly calibrated virial BH mass re- 
lations for Narrow Line Seyfert 1 galaxies and we have shown 
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that, after taking into account radiation pressure, those galax- 
ies seem to have BH masses similar to that of other broad line 
AGNs and and which follow the same MBH-CeXsph relations 
as normal galaxies. 

The small BH masses previously found in NLS 1 can be at- 
tributed to the neglect of radiation pressure in objects radiat- 
ing close to their Eddington limit. 

Overall, the analysis presented in this paper clearly indi- 
cates that radiation forces arising from the deposition of mo- 
mentum by ionizing photon constitute an important physi- 
cal effect which must be taken into account when computing 



virial BH mass estimates. 
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